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Abstract: The normalization of the gluon condensate and of renormalon-related 
power corrections in QCD is computed under the assumption that their "pertur- 
bative" part dominates over any eventual extra contribution from the non-trivial 
vacuum. The calculation is performed in the infrared finite coupling framework, as- 
suming an infrared fixed point is present in the perturbative coupling down to low 
values of Nf. The freezing perturbative coupling is reconstructed using a Banks-Zaks 
expansion approach. Parameter-free predictions of the low energy moments of the 
coupling, which determine the process-independent part of the power corrections, are 
^ \ obtained for a number of choices of the running coupling. 
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1. Introduction 

Recent years have witnessed a rapid development of the area covering the border- 
line between "perturbative" and "non-perturbative" physics in QCD. In particular, 
perturbative ideas have been pushed increasingly far towards the low-energy fron- 
tier to deal with the phenomenology of power corrections. Still in these advances 
the normalization of power corrections is usually considered as an incalculable "non- 
perturbative" parameter, to be fitted from the data. This has been the situation ever 
since their original introduction U by Shifman, Vainshtein and Zakharov (SVZ). In 
this note, I suggest a possibility to compute these parameters from first principles, 
for the limited class of "renormalon-related" power corrections which include in par- 
ticular the "gluon condensate" . Among the various methods devised to deal with 
these contributions, the "infrared finite coupling" approach [|, [!J stands out as 
a particularly attractive scheme. A framework where this approach can be justified 
has recently been suggested 0, and the possibility of a calculation of power cor- 
rections from perturbative input has been pointed out. The aim of this paper is to 
implement the latter suggestion. After a brief review (section 2) of the infrared (IR) 
finite coupling approach and of the proposal in |J, a method to construct the IR 
finite coupling from the Banks-Zaks expansion is described in section 3. The results 
for the power corrections are given in section 4, and further discussed in section 5 
which also contains the conclusions. 
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2. A framework for the IR finite coupling approach to power 
corrections 



In this approach the power corrections are parametrized in term of low energy mo- 
ments of a "non-perturbative" coupling a = a + 5a assumed to be IR finite, where 
a (a = ^ ) is the "perturbative part", and 5a a "non-perturbative" modification 
needed to make a IR finite. Although the approach can also deal with Minkowskian 
quantities, consider as a simple example the case of an Euclidean observable D(Q 2 ) 
in the "single dressed gluon exchange approximation" 

f°° dk 2 




where Q(k 2 /Q 2 ), the "loop momentum distribution function" 0, is known || from 
the relevant single dressed gluon diagrams. Introducing an IR cut-off /ij to separate 
long and short distances, the right hand side is approximated at large Q 2 by 



(2.2) 



with A n (/i 7 ) = tf 1 f£ a(k 2 ) k 2n , where I assumed that <&{k 2 /Q 2 ) ~ c n (k 2 /Q 2 ) n at 
low k 2 and the contribution of the 5a piece has been neglected above /xj. These steps 
can now be justified with the following two crucial assumptions: 

i) The perturbative part a of the coupling reaches a non-trivial IR fixed point 
at low scales and is IR finite by itself, without the need for an hypothetical 5a 
contribution. This statement is likely to be correct for Nf slightly below 16.5 where 



the perturbative coupling has a Banks-Zaks fixed point |L0], [IT], [12] beyond one-loop, 
and I assume it is still true down to Nf = 0. This assumption is supported |13], [14|, [L5| 
by the behavior of the Banks-Zaks expansion for some QCD effective charges. 

Actually, the previous statement must be correct within a range NJ < Nf < 16.5 
which defines the "conformal window" where the perturbative coupling is IR finite 
and causal [|1^, ||. Within the conformal window, there is by definition no 5a term, 
and we have 

D(Q 2 ) = D^iQ 2 ) (2.3) 

with 

^(^ff <* 2 >*(f) < 2 - 4 > 

At large Q 2 , one obtains as in eq. fl2.2|) 
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with the normalization of the power correction 



rrt dk 2 

C PT ^i) = c n j Q — a(k 2 )k 2n (2.6) 

given by a low energy moment of the perturbative coupling. Since the latter is no more 
causal below NJ (even if it still IR finite there), eq. (|2.3| ) cannot be correct anymore 
for Nf < NJ, where the "conformal window amplitude" D-p7p(Q 2 ) is expected to have 
unphysical Landau singularities in the (complex) Q 2 plane. We must therefore have 

D(Q 2 ) = D TT {Q 2 ) + D WP {Q 2 ) (2.7) 

where 1 the "genuine non-perturbative piece" Dj^p{Q 2 ) cancels the Landau singulari- 
ties present in DjPp(Q 2 ). In the standard IR finite coupling approach this piece would 
correspond to the contribution of the Sa part of the coupling in eq.( [2.1|) . Since the 
existence of such a term is quite hypothetical, I shall not assume that the Djfp(Q 2 ) 
piece is related to a (universal) non-perturbative QCD coupling. Still at large Q 2 this 
piece may contribute a "non-perturbative component" Cnp to the 0(l/Q 2n ) power 
correction 



Dn?{Q 2 ) - (2-8) 



so that below NJ we have 



D(Q 2 ) ~ D^Q 2 ) + ^ (2.9) 



hence 



(2.10) 



with 



C(m) = Cp T {h i ) + C np (2.11) 

ii) The second crucial assumption I shall make is that the "non-perturbative" 
component Cnp can in fact be neglected (for not too small fij) in eq. ([2.11|) . This 
assumption, which actually takes the exact counterpart of the SVZ hypothesis [|TJ 
that the "genuine non-perturbative piece" Cnp dominates over the "perturbative" 
fluctuations, can be justified in a number of ways. One is to observe that Cnp, 
which vanishes identically for NJ < Nf < 16.5 within the conformal window, may still 
be small for Nf < NJ below the conformal window, provided Nf is close enough to 

1 If Dpy(Q 2 ) is interpreted as the analytic continuation in Nf of the full conformal window 



amplitude, the decomposition eq.(2.7) is general M, n] and valid beyond the single dressed gluon 



approximation of eq.(2.1) 
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Nj. In fL6| , ||, it was found that in fact 4 < N* < 6, which makes it at least plausible 
the neglect of Cnp at the "real life" QCD value Nf = 3. Another (more drastic) 
possibility is that the power corrections in D-p^iQ 2 ) and Djj-p(Q 2 ) do not match 
(even though at low Q 2 the two components cancel their mutual Landau singularities 
below the conformal window), i.e. that the power corrections are either entirely 
"perturbative" and contribute only to Dpr(Q 2 ), or entirely "non-perturbative" and 
contribute only to Djjp(Q 2 ). This would mean that C^p = even below NJ, and 
only the Cpx component is present, for those condensates (like the gluon condensate) 
which do not vanish within the conformal window, whereas Cnp ^ below NJ only 
for those condensates (like the quark condensate) which vanish identically within 
the conformal window, and therefore have no Cpt component. In such a case, the 
neglect of C^p would be justified at all N/s for the "conformal window type" of 
power corrections. Anyway, the working hypothesis in the following shall be that 
one can compute the bulk of the latter type of power corrections from eq.( [2.6| ) alone. 
In this way, the IR finite coupling approach not only finds a natural framework, but 
its predictiveness is enhanced since there is not any more any "non-perturbative" 
free parameter and the normalization of power corrections can be computed, as we 
demonstrate in the next section (in this sense the approach goes beyond the operator 
product expansion even when applied to Euclidean quantities). 

3. Reconstructing the IR finite perturbative coupling: a Banks- 
Zaks expansion approach 

Even though the perturbative coupling appears to have an IR fixed point for large 
enough Nf beyond two- loop, this is not always manifest when one decreases Nf. For 
instance, the Banks-Zaks fixed point at two-loop relies on having ft < 0, which is 
not realized for Nf < 8. Then one might rescue the fixed point with a negative 
three-loop term, but even this feature is usually lost at Nf = 3. On the other hand, 
as mentioned in section 2, the Banks-Zaks expansion does signal in a number of cases 
the persistence of the fixed point even down to Nf = 2. This observation suggests the 
following strategy: try to reconstruct the IR finite coupling, and eventually supply 
the missing higher order terms in the beta function (5(a), given the IR fixed point 
Banks-Zaks expansion. This is an expansion in powers of the distance 16.5 — Nf 
from the top of the conformal window, which is proportional to ft. The solution 
a* = a*(e) of the equation 

[3(a) = fta, e) = -ft a 2 - ft a 3 - ft a 4 - ft a 5 + ... = (3.1) 
in the limit ft — ► 0, with ft [i > 1) finite is obtained as a power series 

a* = a *(e) = e + 5 1 e 2 + 5 2 e 3 + ... (3.2) 
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where 



x a #2,0 
0i = Pi,i - -5— 

Pl,0 

82 = #1 + 9i Si + P 2 ,i ~ 92 (3.3) 

The expansion parameter is e = ^-(16.5 — Nf) = —-j^- The (3 it j, which are 
Af/-independent (but scheme dependent for i > 1), are defined by /3i = /?i,o + Pi,i Po 
(pi,o = ~if , = x)"^2 = #2,0 + #2,1 Po + p2,2 ^0 (I assume f3 2 is at most quadratic 
in Nf, hence in f3 ) and g\, g 2 are given in eq.( |3.6| ). Given the knowledge of the 3-loop 



beta function in (e.g.) the MS scheme, 02,o can be obtained |T7[ from a one-loop 



calculation of a (see eq.( |4.2|) ). I shall also use the related expansion for the critical 
exponent 



7 = 7(e) = fV,e) (3.4) 
= -(3 lfi e 2 (l + gi e + g 2 e 2 + ...) (3.5) 

where 



91 = Pi,i 

2 . P3,0 ( $2,0 \ 

92 = 9l + A^-fej 

The giS are scheme independent |12| , and |13| g 2 = —8.89. 



(3.6) 



The method relies on the differential equation [12| for a*(e 



(3.7) 



where via, e) is the Nf dependent part of the beta function, after splitting off the 
e = (i.e. Nf = 16.5) piece /3(a,0) 



/3(a, e) = (3{a, 0) — e v(a, e) 
Its expansion in powers of a is 



(3. 



u(a, e) = -p lfi a 2 [l + (3 hl a + (/3 2il + /3 2j2 (3 )a 2 + ...] (3.9) 

Eq. (|3.7|) follows by taking the total derivative with respect to e of the relation 
/3(a*, e) = which defines the fixed point a*(e) 



d(3 da* dp 







(3.10) 
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and using eq.Q) and § (a, e) = -|[e i/(a,e)] (eq.Q). 

It is convenient to introduce the function e*(a), which is the inverse of the Banks- 
Zaks function a*(e): for given a, e*(a) is the value of e (i.e. of iV/) where (3(a, e) = 0. 
The knowledge of e*(a) and of z/(a, e) determine /3(a, 0), hence the full beta function. 
Indeed using eq. fl3.8Q the condition (3[a,e*{a)] = becomes 



P(a,0) = e*(a) u[a,e*(a)} (3.11) 

Hence 



P(a, e) = e*(a) u[a, e*(a)} - e v(a, e) (3.12) 
In term of e*(a) eq. (|3.7| ) reads 

fl[ e „(o,e)]|«. = ^ (3.13) 

da 

Eq. (|3.13Q gives a constraint on z/(a, e) given the Banks-Zaks functions 7(e) and a*(e). 
This constraint is not sufficient to determine u(a, e) (and the beta function) without 
further assumptions. In the following I shall assume that z/(a, e) = z/ (a) is indepen- 
dent of e, i.e. that the beta function coefficients are at most linear in Nf (or /3q): 
this amounts to an approximation, in the spirit of the Banks-Zaks approach, where 
one keeps only the leading e = term in an expansion of u(a, e) in powers of e (in 
particular, one neglects the $2,2 Po term in eq. (|3.9Q ). Then eq.( p,13| ) gives 

Ma) = (3.14) 

da 

and from eq.( p,12|) one gets 



P(a,e) = [f{a)-e] (3.15) 

da 

Using the Banks-Zaks expansions of the fixed point a*(e) and of the critical exponent 
7(e) truncated to a given order as input, eq.( |3.15[) yields a corresponding "improved" 
approximation to the beta function, which displays a built-in fixed point at a = a*(e). 

In this approach, the leading order (LO) approximation thus gives e*(a) = a and 
7[e*(a)] = — /5i,o a<2 - The next-to-leading order (NLO) approximation uses the NLO 
Banks-Zaks expansions of the fixed point and of the critical exponent: e*(a) is then 
obtained by inverting eq. (|3.2j ) (with 62 = 0), i.e. solving for e* in a = e* + 5\ e* 2 and 
reporting in eq. (|3.15Q , with j(e*) = — /3 1)0 e* 2 (l + gi e*). The next-to-next-to-leading 
order (NNLO) approximation uses the NNLO Banks-Zaks expansions of the fixed 
point (which requires the knowledge of /?2,i) and of the critical exponent (eq.(|3.2|) 
and (B.5I)), etc... 
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The approximation can be further systematically improved by including the 
knowledge of the known iVf-dependent terms in the beta function. For instance, 
if the three-loop fa coefficient is known, one can include the knowledge 2 of the term 
quadratic in f3 Q in /3 2 with the ansatz 

u(a, e) = v {a) + (3 2 fi f3 2 ,2 a 4 e (3.16) 

where vo(a) = v(a, 0) is independent of e (the knowledge of @2,o and /?2,i is contained 
in the NLO and NNLO terms in the Banks-Zaks expansion of a*, as mentioned 
above). Eq. (|3.13| ) then fixes ^o(a) from 

vo{a) + 2 ft iQ f3 2i2 a 4 e*(a) = (3.17) 

da 

which yields u(a, e), hence from eq. ( |3.12|) 

/3(a, 6) = [e* (a) - e] - (3 2 1>0 /3 2 , 2 a 4 [e*(a) - e] 2 (3.18) 



4. Results 



The "coupling" appearing in eq.fl2.4p should be viewed as a physical, gauge-independent 
quantity, just as the observable D(Q 2 ) to which it is directly related. In the IR finite 
coupling approach, it is also assumed to be universal, i.e. the same for all observ- 
ables. The existence of such an object is still speculative. It is attractive to identify 
this coupling to the "skeleton coupling" [0, |T8| associated to a (yet hypothetical) 
"QCD skeleton expansion". A promising approach in this direction is provided by 



the "pinch technique" construction [19], BO]. The pinch coupling is presently known 



only at one- loop, where it is related to the MS coupling by 



<k 2 ) = aW s(fi 2 ) + [-/3 (log(A: 2 //i 2 ) - 5/3) + d t ] af^(/i 2 ) + ... (4.1) 

with di\pinch = 1. An alternative suggestion [||] is to use the "gluon bremsstrahlung 
coupling" ||21|1 , also known to the one- loop level eq. (|4.1| ) with di\brems = 1 — 7r 2 /4. 
Since the full three-loop beta function coefficient (hence 82) is not yet known for 
these two couplings, I shall apply the method of section 3 in the NLO approximation 
described there. Actually, since the Banks-Zaks expansion of the critical exponent 
is known [TjJ up to NNLO (eq.( |3.5| )), and may be reliable || even down to Nf = 3, 
I shall go half-way towards the NNLO approximation, and use eq. fl3.5p in eq. (|3.15| ), 
while still using eq.fl3.2p (with 82 = 0) to fix e*(a). The input scheme dependent 
numerical values following from eq. fl4.1| ) and the relation [ jTTH 

2 Actually, given that < a < a* = 0(e), this term is effectively of the same order as the f3^,i a 5 
term in v(a,e) (eq.(3.9)), and should be taken as input only together with the latter, i.e. at the 
NNNLO level. 



7 



are @ &&\pinch = 2.61 and ^\brems = 2.61 + ^ = 5.08. Hence 5i|pmc/i = 2.14 
while 5\\brems = —0.33 (a smaller correction!). It follows from eq. fl3.2|) that at 
Nf = 3 the IR fixed point a* = 0.299 for the gluon bremsstrahlung coupling, smaller 
then the corresponding value a* = 0.578 for the pinch coupling which is subject to 
rather large uncertainties. 

As a third alternative, I would like to suggest the "universal coupling" introduced 
in fl2"l , because of its simplicity. It is defined 3 by the condition 5i = for all z's, i.e. 
a*(e) = e, and therefore its beta function can be expressed entirely in term of the 
critical exponent 



/3(a, e) = (a — e)j(a) 

= -(A, a 2 + (3 1>0 a 3 )(l + gia + g 2 a 2 + ...) (4.3) 

At Nf = 3 the IR fixed point is a* = e = 0.336. The scale is fixed knowing that 
j^\universal = /3i 5 i (from 5i = 0), which determines (eq.( |4.2|) ) di\universal = 
—1.137, and the natural assumption that the term proportional to /3q in eq. (|4.1|) is 
the same. 

The bremsstrahlung coupling beta function at Nf = 3 is shown in Fig.l. Note the 
negative ultraviolet fixed point at a ~ —0.17. It corresponds to a zero of the critical 
exponent eq. (|3.5| ) at e ~ —0.15, and is a necessary condition for the scenario in || 
to determine the bottom Nf of the conformal window from the condition 7(e) = 1, 
which yields Nf ~ 4 if one uses eq.( |3.5|) . The resulting running coupling is shown 
in Fig. 2, where I used af rs (Mz) = 0.117 as input (eq.(£OJ) yields the corresponding 
input value of the bremsstrahlung coupling) . 

It is then straightforward to compute the first few low energy moments of the 
coupling 

a 2n -iM = J Q n — \ -^\ a(k 2 ) (4.4) 
In term of the beta function they are given by 

a-Tn-iM = n ~7T/ r a exp[n$(a,a 7 )] (4.5) 

J a* 

where 

3 I assume linear Nf dependence. Otherwise there is the more general solution /3(a, e) = 7(0) (a — 
e)+7i(a)(a-e) 2 + .... 
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Figure 1: The bremsstrahlung coupling beta function (Nf = 3). 




GeV 



Figure 2: The bremsstrahlung coupling (Nf = 3, a MS (M z ) = 0.117) 



^ a ^Lwr Xo6 K?) (46) 

is the solution of the renormalization group equation with aj = a(fij). Taking 
fij = 2GeV, one gets the results in Table 1 if af ts (Mz) = 0.117 and those in Table 
2 if af s (M z ) = 0.120. Note the sensitivity to the high energy input value of a s . 

These results are subjected to theoretical uncertainties, stemming from the mag- 
nitude of the IR values of the coupling which should induce seizable higher order cor- 
rections. The convergence of the Banks-Zaks expansion is bad in the pinch coupling 
case (which has a large IR value), and the knowledge of the 3- loop beta function 
coefficient and of 5 2 is essential for a more reliable prediction. The situation looks 
better for the "universal coupling" and the bremsstrahlung coupling. To assess the 
convergence of the expansion, the results for the moments in the NLO approximation 
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where one uses only the first two terms in the Banks-Zaks expansion of 7 (eq.( |3.5| )) 
are quoted within parenthesis in the tables. 

The n = moment gives the process-independent part of the normalization of 
the 1/Q power corrections. If one uses the gluon bremsstrahlung ansatz for the 
coupling, the predicted value is in qualitative agreement with the experimentally 
determined f| one (a ~ 0.14 — 0.17), although it should be remembered that the 
latter depends on the way the "perturbative part" of the amplitude (the piece above 
Hi in eq.( [2.10| )) is handled, as well as upon extra assumptions in the case of non- 
inclusive Minkowskian observables. The n = 3 moment gives the normalization of 
the "gluon condensate" 

< —G 2 > w = i^/i} (4.7) 

Note the definition used here involves an arbitrary IR cut-off /i/, as necessary in 
the case of renormalon-related power corrections. If one wants to compare 4 to the 
effective phenomenological SVZ definition PJ, one can just compute the integral in 
eq. (|2.4j ) (which does not depend on /i/) for any given Euclidean observable where the 
gluon condensate gives the leading power correction, and fit the result with the SVZ 
ansatz. For instance, for the Adler D function 

D(Q 2 )^a(Q 2 ) + 2 -f±<^G 2 > (4.8) 

where the ^zd factor is the leading order coefficient function. Similarly, the SVZ 
condensate < ^-G 2 > could be defined from the basic observable 03(^7) (eq. (|4.4j ) 
with n=2), where the IR cut-off fii now plays the role of the high energy scale Q, by 

a 3 (/ij) ~ a(fxj) + ^-— A < —G 2 > (4.9) 

For fu = 2GeV, eq.flOD yields < ^G 2 >~ 0.05 GeV 4 for crp(M z ) = 0.117, which 
looks reasonable compared to the standard SVZ value. However, this comparison 
is actually devoid of significance due to the following intriguing fact: varying the 
scale /if in eq.(|4.9j), one finds the discrepancy 03(^7) — a(/if) between 0,3(^1) and 
its lowest order perturbative approximation a(fij) decreases much slower then the 
inverse fourth power of hj\ A similar result is obtained if one uses eq. ([4.8|) (with 
the "loop momentum distribution function" Q(k 2 /Q 2 ) taken from 0). Since the 
(principal value regulated) Borel sum of the perturbative series associated to the 
observables a 3 (/i/) (or D(Q 2 )) are known |22|, |23|] to differ from the exact values by 
just such an 0(1/ /ij) (resp. 0{1/Q A )) correction, one is bound to conclude that the 
naive treatment of approximating the Borel sum by its leading order term does not 

4 I am indebted to Al. Mueller for raising the question. 
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work 5 here. This is another point of discrepancy with the standard SVZ procedure, 
on top of the assumption that the "perturbative part" of the condensate dominates. 



5. Discussion and conclusions 



The essential assumption in the present approach is that the perturbative beta func- 
tion has an IR fixed point at least down to Nf = 3. This is partly implemented 
by constructing beta functions with negative three-loop coefficients: at NNLO the 
method of section 3 yields /3 2 = p 2 ,o + /?2,i Po where both /? 2 ,o an d /?2,i turn out 
negative (see footnote 4) for the considered couplings. Actually, essentially the same 
results can be obtained (at least for the bremsstrahlung coupling) in a simpler way, 
which makes it transparent the reason for the existence of the IR fixed point. In- 
deed, consider the 4- loop beta function eq.(|3.1p, and observe that in the IR region 
the usual power counting should be modified: namely, given that a is 0((3q) there, 
to O (a 6 ) accuracy one should drop the 0((3q) term in /3 2 , and keep only the leading 
0((Pq)°) term in /3 3 , i.e. use the effective 4-loop beta function (in accordance with 
the remark in footnote 1) 



P eff {a, e) = -fa a 2 - fa a 3 - (/3 2 , + Po)a 4 - /3 3 ,o a 5 + 0(a 6 ) (5.1) 

(in the ultraviolet region, this beta function has of course only the 0{a i ) accuracy 
of the 2- loop beta function). In the case of the bremsstrahlung coupling, the results 
obtained using the 4-loop /3 e ff turn out to be very close to those of section 4 in 
the NNLO approximation. For instance /3 e // has an IR fixed point at a* = 0.294 if 
Nf = 3 (I used /3 3) o = 37.76 from eq. (|3.6p ), and one gets: ao = 0.201, a\ = 0.177 
and_a 3 = 0.156 if af^(M z ) = 0.117, and a = 0.210, a l = 0.189 and a 3 = 0.168 if 
a ! f s (Mz) = 0.120. Similarly in NLO one should use a 3-loop /3 e //(a, e) = — (3q a 2 — 
Pi a 3 — /3 2 ,o ° 4 + C(a 5 ), and in LO a 2-loop P e ff(a, e) = — (3 a 2 — /3 lj0 a 3 + 0(a 4 ). 
The presence of an IR fixed point in j3 e ff down to low values of Nf seems to be a 
general phenomenon, at least up to NLO. This is obvious in LO, since (3\q is scheme 
independent, but less so in NLO where /3 2) o is scheme dependent. Nevertheless it 
turns out that /3 2j0 is negative for all known physical effective charges []T8"|, as well 



as for the three couplings quoted above. Consequently, there may be a positive zero 
in the 3-loop (3 e ff correctly signalling an IR fixed point, even if the standard 3-loop 
beta function has no positive zero with all its coefficients of the same sign. 

At NNLO, the presence of an IR fixed point in the 4-loop /3 e // may be jeopardized 
by large positive values of /?2,i and (or) /3 3i0 . Actually, /3 2 ,i turns out to be negative 
for all known 6 effective charges (except the one ("ay") defined by the static QCD 



5 Similar results are obtained if one uses || the BLM scale Eq] in a. 

6 For the pinch coupling and the bremsstrahlung coupling, ^2,1 nas been "predicted" from the 
assumption that 82 — 0, which yields (eq.(|3.3|)) /?2,i = —23.6 for the pinch coupling and /?2.i = 
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potential, where it is positive |H| but small enough not to destabilize the fixed 



point). The real problem comes from the 4-loop coefficient /^q, which is positive for 
all known effective charges (except again ay, where it is negative and tiny). In the 
case of the pinch coupling, it turns out in fact too large (one gets /3 3i0 = 164.7 from 
eq.( |3.6| )) for the 4-loop (3 e ff to have an IR fixed point if Nf < 13. Similarly, in the 
case of the Adler D-function effective charge where fi^o = 127, the 4-loop j3 e ff does 
not have an IR fixed point if Nf < 11. For all other effective charges however the 
4-loop (3 e ff does exhibit an IR fixed point down to Nf = 0! However, in those cases of 
large positive /3^o (which is a consequence of a small 02,0, see eq. (|3.6j )), the method 
of section 3 provides an effective resummation of the relevant higher order terms, 
obtained under the assumption the Banks-Zaks expansions of the IR fixed point and 
of the critical exponent do converge: all known effective charges then appear 7 after 
resummation to have an IR fixed point down to Nf = (although the convergence 
of the fixed point Banks-Zaks expansion becomes problematic already at Nf = 3 for 
some of them, such as the pinch coupling). 

The suggestion of perturbative freezing of the coupling at low Nf was first made 
in [T3]. There is however an essential difference with the present proposal: it is not 



suggested here that the perturbative IR fixed point has anything to do with the low 
energy behavior of the full QCD amplitudes below the conformal window, which is 



entirely non-perturbative. For instance, as observed in [g5j spontaneous chiral sym- 
metry breaking considerations at large N c imply the Adler D-function must vanish 
at zero momentum, which is inconsistent with the positive value expected from per- 
turbative freezing. What is suggested instead is that the perturbative freezing at 
low Nf is relevant to determine the normalization of renormalon- related condensates 
and power corrections which appear in the short distance expansion of amplitudes. 
This amounts to the recognition that objects like the "gluon condensate", at the dif- 
ference of the quark condensate, are of a basically "perturbative" nature, and thus 
unrelated to "genuine" non-perturbative properties of the vacuum such as chiral sym- 
metry breaking or confinement. The notion of a "conformal window" is an essential 
part of the present proposal: only those power corrections which are already present 
within the conformal window are amenable to a perturbative treatment, and below 
the conformal window there are other really "non-perturbative" contributions which 



—7.43 for the bremsstrahlung coupling. This assumption turns out to yield rather good results 
in the case of the effective charges associated to the Adler D-function and the polarized (gi) and 
non-polarized (Fi) Bjorken sum rules, for which the "predicted" values are respectively ,#2.1 = 
— 16.17,-9.98,-6.97 compared to the exact values (corrected for some numerical inaccuracies in 
Jig] ) fa.i = —15.94,-11.19,-6.81. The partial reason for this success are the large cancellations 
between gi (the "scheme independent" contribution to 62 in eq.( |3.3| )) and the "scheme dependent" 
contribution which involves S\ and 02,1 ■ 

7 This is even true for the effective charge associated |24[ to Higgs decay. In this case however one 
gets a large fixed point value a* = 0(1), and convergence of the Banks-Zaks expansion is doubtful 
for N f < 4. 
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are crucial to determine the true low energy properties of QCD. Moreover, it was 
shown in || that the assumption that the perturbative IR fixed point persists below 
the bottom N* of the conformal window leads to the condition 7(e) = 1 to determine 
Nj. It is interesting that this condition gives NJ ~ 4, rather close to the "real life" 
QCD value Nf = 3, which might give an alternative justification to the suggested 
calculation procedure based on the "anti-SVZ" hypothesis that the "perturbative" 
piece of the condensates actually dominates over the "non-perturbative" fluctuations. 
Note that the opposite SVZ assumption of dominance of the non-perturbative piece 
has been questioned previously in the literature (see e.g. @). Furthermore, even 
if the present assumption turns out to be invalid, the results of this paper are still 
useful to extract from experiment the "truly non-perturbative" part Cnp, which is 
given a completely unambiguous definition through eq. ( |2.9| ). 

The typical example of a "perturbative" conformal window amplitude is the "sin- 
gle dressed gluon" integral of eq. ( j2.4| ) , where the running coupling inside the integral 
is IR finite, and calculable from perturbative input through an (eventually resummed) 
Banks-Zaks expansion. It is implicitly assumed that this particular coupling is free of 
IR renormalons and can be unambiguously determined from its perturbative series 
(say, by Borel summation). The corresponding Banks-Zaks series should then be 
also Borel summable. Note also that the integral eq. Q2.4|) is free of any renormalon 
ambiguity, although renormalons are present in the corresponding perturbative se- 
ries, but is still expected to be affected below the conformal window by unphysical 
Landau singularities in the complex Q 2 plane. Such an amplitude |23[ represents 
a natural form of a generalized perturbation theory, which gives the background on 
top of which genuine non-perturbative contributions may take place below N*. The 
calculation of the "perturbative condensate", although using only perturbative in- 
formation, goes beyond a mere renormalon estimate, since there is usually a part 
[f2~6| , H of the low momentum contribution of the perturbative coupling which is not 
determined only [^6, 27] by renormalons. The assumption that renormalon- related 



power corrections are "perturbative" in the above sense also gives a straightforward 
justification to the IR finite coupling approach to power corrections, and leaves no 
arbitrary free parameter (except of course the overall QCD scale) to be fixed from 
experiment. The main conceptual problem in this framework remains to find the 
identity of the (hopefully unique) perturbative IR finite QCD coupling which deter- 
mines the power corrections, and derive the systematic form to all orders of the (yet 
to be constructed) generalized perturbation theory, perhaps |18| along the lines of 
a QCD "skeleton expansion" : this is however a problem of a basically perturbative 
nature. 

Acknowledgments 

I thank Yu.L. Dokshitzer and A.H. Mueller for stimulating discussions. I also wish 



13 



to thank the referee for constructive criticism. 



References 

[1] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B 147 (1979) 385. 
[2] See e.g. M. Beneke, Phys. Rept. B 317 (1999) 1 |hep-ph/9807443|. 



[3] Yu.L. Dokshitzer and B.R. Webber, Phys. Lett. B 352 (1995) 451 ||hep-ph/95042ig| 1 . 



[4] Yu.L. Dokshitzer, G. Marchesini and B.R. Webber, Nucl. Phys. B 469 (1996) 93 ]Kcp~ 
ph/951233^. 



[5] For a review and further references, see Yu.L. Dokshitzer, in International Conference 
"Frontiers of Matter", Blois, France, June 1999 [ |hep-ph/99lT299| . 



[6] G. Grunberg, Conformal window and Landau singularities, hep-ph/0104098 . 
[7] G. Grunberg, J 'HEP 11 (1998) 006 [|hep-ph/9807494 1 . 



E. Gardi and G. Grunberg, JEEP 11 (1999) 016 [|hep-ph/9908458|| . 



[9] M. Neubert, Phys. Rev. D 51 (1995) 5924 [ |hep-ph/ 94 12265 

[io; 



[11 

[12 

[13! 

[14 
[15 

[16! 
[17! 
[18! 

[19 

[20! 
[21 



T. Banks and A. Zaks, Nucl. Phys. B 196 (1982) 189. 

A.R. White, Phys. Rev. D 29 (1984) 1435 ; Int. J. Mod. Phys. A 8 (1993) 4755 
[|hep-th/9303053 |. 

G. Grunberg, Phys. Rev. D 46 (1992) 2228 . 



P.M. Stevenson, Phys. Lett. B 331 (1994) 187 ||hep-ph/9402276|| ; S.A. Caveny and 
P.M. Stevenson, |hep-ph/9705"3l9. 



E. Gardi and M. Karliner, Nucl. Phys. B 529 (1998) 383 [frep-ph/9802218 |. 
E. Gardi and G. Grunberg, JEEP 03 (1999) 024 ||hep-th/98101~92 |. 



G. Grunberg, Fixing the conformal window in QCD, hep-ph /0009272 . 
G. Grunberg, Phys. Lett. B 135 (1984) 455. 

S.J. Brodsky, E. Gardi, G. Grunberg and J. Rathsman, Phys. Rev. D 63 (2001) 094017 
|iep-ph/0002065 |. 

N.J. Watson, Nucl. Phys. B 494 (1997) 388 [[hep-ph/960638l|| ; Nucl. Phys. B 552 
(1999) 461 [|hep-ph/9812"202|; |hep-ph/9912303|. 



J. Papavassiliou, Phys. Rev. D 62 (2000) 045006 ||hep-ph/9912338|| ; Phys. Rev. Lett. 
84 (2000) 2782 jhep-ph/9912336|| . 

S. Catani, G. Marchesini and B. R. Webber, Nucl. Phys. B 349 (1991) 635. 



14 



[22] G. Grunberg, Phys. Lett. B 372 (1996) 121 [|hep-ph/9512203;l ; in DIS96 fhep- 



ph/9608375|] . 



[23] Yu.L. Dokshitzer and N.G. Uraltsev, Phys. Lett. B 380 (1996) 141 |hep-ph/9512407 |. 

[24] J. A. M. Vermaseren, S.A. Larin and T. van Ritbergen Phys. Lett. B 405 (1997) 327 
Jhep-ph/9703284 |. 



[25] S. Peris and E. de Rafael, Nucl. Phys. B 500 (1997) 325 |hep-ph/97014T8 ]. 
[26] G. Grunberg, Phys. Lett. B 325 (1994) 441. 
[27] G. Grunberg, Phys. Lett. B 349 (1995) 469. 



15 





bremsstrahlung 


universal 


pinch 


a 


0.207 (0.222) 


0.225 (0.246) 


0.330 (0.366) 


ai 


0.176 (0.198) 


0.187 (0.212) 


0.256 (0.300) 


a 3 


0.155 (0.173) 


0.163 (0.183) 


0.210 (0.243) 



Table 1: Moments for af s (M z ) = 0.117. 





bremsstrahlung 


universal 


pinch 


a 


0.217 (0.232) 


0.237 (0.259) 


0.353 (0.390) 


a\ 


0.188 (0.213) 


0.202 (0.230) 


0.281 (0.332) 


a 3 


0.167 (0.189) 


0.176 (0.201) 


0.233 (0.273) 



Table 2: Moments for af S (M Z ) = 0.120. 
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